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The mathematical relationships that describe an amplitude-modulated 
data system are developed in a form suitable for programming on a high- 
speed digital computer. These equations contain expressions that specify in 
general terms the transmission-frequency characteristics of a transmission 
medium. A data signal composed of a train of raised-cosine shaped pulses 
is generated in the stimulating process. The simultation provides a means 
for computing the resulting response of systems to pulse trains. The per- 
formance of a double-sideband AM data system is evaluated from meas- 
urements of the maximum vertical opening, or aperture, of the eye pattern 
formed by the received signal. This aperture is related to the system per- 
formance in terms of signal-to-noise ratio and error rate of the system. A 
verification of this technique is made by simultaing the conditions of an 
experimental laboratory data system on the computer and comparing com- 
puted and measured performance. 

I. INTRODUCTION 

The process of transmitting digital information and the problems re- 
lating to this type of communication have received much attention in 
recent years. Practical means of communication based on this type of 
information transmittal have been in use for a long time in the form of 
telegraph and teletype systems, but for both these systems the speed of 
operation is relatively slow, since the information is supplied either by a 
manual operator or by a mechanical device of some sort. Recently, how- 
ever, many new sources of digital information have arisen for which the 
information rates may range from telegraph speed up to several orders 
of magnitude greater than that speed. This has created a need for new 
high-speed data transmission systems. To this need for higher speed has 
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been added the requirement of greater accuracy — amounting to, in 
some cases, essentially error-free transmission. For a given bandwidth 
both high speeds and accurate transmission are often difficult to realize. 
For this reason, a great deal of effort has been spent in devising means 
for utilizing a channel as effectively as possible in order that both speed 
and accuracy may be maximized. 

This effort may be divided into two areas. The first is concerned with 
what may be called the terminal problem, in which effort is directed to- 
ward developing data system terminal equipment that makes best use of 
the transmission channel. This has resulted in a number of competing 
schemes for data transmission based on different modulation methods, 
such as double and vestigial sideband amplitude modulation, frequency 
modulation and phase modulation. Also, a number of different detection 
methods are available for these various systems. 

The second area of effort may be referred to as the transmission 
problem. Here, work is directed toward determining what influence the 
various transmission factors such as the transmission-frequency char- 
acteristics and noise have on the transmitted data signal, and what 
conditions should be maintained for a particular data system that is to 
operate at a given speed and accuracy. 

In order to answer questions that arise in these two areas, it is neces- 
sary that a method for evaluating data system performance be devised. 
Good performance is associated with both high speed and accurate 
transmission, but quality of performance is not necessarily preserved by 
a direct interchange of these two quantities. It is of little value to send 
information at very high speeds with very low accuracy. The optimum 
performance is obtained, therefore, by transmitting data over a given 
channel at the highest rate possible, consistent with a given error cri- 
terion as dictated by the system requirements. 

Because of the many factors that contribute to the over-all perform- 
ance, the problem of evaluating performance is a difficult one. The most 
direct method and the one that has been most widely used is that of 
measuring performance on a data system under actual operating condi- 
tions. Some results from such experimental tests have been given 1 and 
made general, at least to some extent. The experimental results, however, 
are of necessity limited in their generality, and the task of obtaining 
results by this method is expensive and time-consuming. In addition, 
for design purposes it is desirable to know how a data system will per- 
form before it is built as well as after. 

Another approach to evaluating performance is by means of a theo- 
retical analysis of the problems involved. A considerable amount of work 
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has been done in this area. 23 - 45 While this work is of great value, the 
complexity of the problems make a general solution untractable, except 
under conditions of rather restricted and idealized assumptions. For 
conditions of practical significance, however, a system analysis using 
mathematical models is useful, provided the task of a numerical computa- 
tion can be carried out. 

This paper is the result of an effort to evaluate the performance of an 
AM data system by methods of numerical analysis and simulation of 
the system on a digital computer. In this work, the general method of 
analysis and simulation is outlined, a criterion for performance is defined 
and some results that have been obtained by the use of this method are 
given. While the problem is analyzed in rather general terms, the par- 
ticular results presented apply to the specific problem of evaluating the 
degradation in performance due to delay distortion. 



II. BASIC DATA SYSTEM 

The basic data system considered here is shown in Fig. 1. The input 
function f(t) is the electrical representation of the random sequence of 
binary information to be transmitted. This signal is made up of a suc- 
cession of identical raised-cosine shaped pulses interspersed with gaps 
or spaces where no pulse occurs, forming a pattern of the type illustrated 
in Fig. 2. In this pattern the occurrence or absence of a pulse is governed 
by the occurrence of a one or a zero in the binary information. These 
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Fig. 1 - Data transmission system. 




110 

BINARY INFORMATION 



Fig. 2 - Binary information and equivalent pulse pattern. 
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binary values correspond in time to the points at which pulses may 
occur in the pattern, and they are separated by a uniform time interval 
called the sampling interval, which is the reciprocal of the bit speed of the 
system. 

The basic pulse used to form the pattern is shown in Fig. 3. It is chosen 
to have a maximum width to , which is equal to two bit intervals at its 
base. Because of this fact, and because of the raised-cosine shape, pulses 
that are adjacent in the pattern of Fig. 2 add together over the interval 
that separates them to give a constant value equal to the maximum value 
of a single pulse. Also, adjacent gaps produce a constant value equal to 
zero. For the case of alternating values of one and zero, transitions are 
produced that are portions of a sinusoid. 

For purposes of analysis it will be useful to consider a single pulse of 
the type used in forming the pulse pattern of Fig. 2. To do this a par- 
ticular pattern is used, as shown in Fig. 3(a), in which /(i) is formed by 
the binary sequence 010 followed by a sequence of zeros and then re- 
peated periodically at intervals of T. This method of specifying the 
signal permits the analysis of the basic pulse to be made on a Fourier 
series basis that is conveniently adapted to numerical computation on 
the digital computer. 




(b) 

Fig. 3 - Basic pulse used to form pattern: (a) raised-cosine pulse; (b) associated 
spectrum. 
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The spectral density function resulting from a Fourier analysis of the 
periodic pulse pattern of Fig. 3(a), is shown in Fig. 3(b), and is ex- 
pressed mathematically as 

C = - (I) 



[■ - (*)T 



We can imagine a perfect data system in which the signal f(t) of Fig. 
2 is transmitted without alteration. This will be considered a reference 
condition for measuring performance. Under this condition the binary 
information can be extracted from the received signal g(t), which is 
identical to /(f), by examining the pulse pattern at the sampling times to 
determine if the signal is greater than or less than some intermediate- 
value called the slice level. (In most cases this level is equal to one half 
the value of the maximum pulse height. Under certain conditions, how- 
ever, it will be advantageous to use other values.) 

For a system with ideal transmission-frequency characteristics, errors 
are produced if a sample is made during the time when a noise disturbance 
is large enough to cause the signal to make an erroneous excursion across 
the slice level. In order to meet system requirements the value of the 
signal-to-noise ratio must be such that the number of errors in a given 
time does not exceed a certain maximum. If the signal-to-noise ratio is 
less than this value the system is considered as failing. If the signal-to- 
noise ratio is such that the error criterion is just met, then a reduction 
in the noise (with the signal held constant) results in an operating margin 
against failure equal to the amount of the noise reduction. If the system 
were operating without noise, the amount of noise necessary to bring it 
to the failing condition would be called the noise margin for the system. 

As the transmission-frequency conditions depart from ideal, the out- 
put g(t) will no longer be identical to the input f(t). This lack of fidelity 
is evidence of system degradation. It will be shown later how this deg- 
radation can be expressed quantitatively as a reduction in the noise 
margin for the system. 

Transmission of the data signal over most facilities requires some form 
of modulation. It will be assumed here that the signal is amplitude- 
modulated by an ideal product modulator, and that an ideal envelope 
detector is used at the receiver. Also, the system will be considered, 
for the present, as being noise-free, and the impairment that the data 
signal suffers in transmission will be due entirely to the frequency char- 
acteristics of the transmission medium. 
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III. TRANSMISSION CHARACTERISTICS 

The basic relationship that describes the transmission-frequency char- 
acteristics of a transmission system is given as 

r(co) = am*-* *, (2) 

in which A(w) is the attenuation and ^(co) the phase characteristic. In 
the work to follow we will make use of the fact that the attenuation 
characteristic has even symmetry A(co) = A (-co) and the phase char- 
acteristic has odd symmetry ^-(co) = — \p( — co) about the zero frequency 
axis. 

To accommodate the use of a Fourier series development of the 
problem, the transmission-frequency characteristics can be expressed in 
terms of their values at discrete radian frequencies corresponding to the 
Fourier series components, rather than on a continuous basis. For this 
purpose, (2) becomes 

T (^\ = A (^A e -m^,T K (3) 



(?) 



T 



The phase curve of Fig. 4(a) has a shape typical of the characteristics 
of many transmission facilities and can be approximated by a sine-wave 
departure from linearity as given by 

\p{u>) = au — b sin (cot + Q), (4) 

where : 

/> is the maximum departure of the phase from linearity, 
r is the reciprocal of the period of the phase curve and 
is the displacement of the phase curve from symmetry about the 

carrier. 

The corresponding envelope delay is shown in Fig. 4(b) and is expressed 

analytically as the derivative of the phase: 

^ = a - br COS (cor + 9). (5) 

ffco 

The linear term in the phase expression results in only a constant 
delay, and the distortion in transmission is therefore produced entirely 
by the sine-wave departure from linearity. The three parameters b, t and 
6 give a complete description of any sinusoidal delay or phase character- 
istic. A wide variety of such curves can be realized by letting these three 
parameters range over values of interest (it will be shown later that, 
with a sinusoidal representation of phase, a very close approximation to 
actual transmission characteristics can be achieved). With this repre- 
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Fig. 4 - Sinusoidal characteristics: (a) phase; (b) delay. 

sentation of the transmission characteristics, the problem of determining 
the performance of a data system in terms of the parameters b, r and 6 
will now be considered. 



IV. ANALYSIS OF AN AM DATA SYSTEM 

In this section,* the analysis of an AM system is made in order to 
obtain the expression for the output function g(t) in terms of the input 
function f(t) and the parameters b, t and 0. For a periodic input pulse 
as shown in Fig. 3(a) the output response will also be periodic. If the 
period T is chosen large enough to allow the transient response of the 
output to be restored essentially to zero there will be no intersymbol 
interference between these periodic pulses. We will then have a condi- 
tion in which both f(t) and g(t) can be represented by a Fourier series. 

The signal that results from an amplitude modulation of f(t) by a 
cosine carrier is given by 



h(t) = f(t) cos w L 



(6) 



* This analysis is based on an unpublished memorandum of R. G. Sogers. It 
includes some slight modification for use in the present study. 
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Expressing f(t) by a complex Fourier series representation and using the 
complex form of the cosine function, we have 

Joint I — jo>a l °° 

h ( t ) = e ±J— 2 C n e jn ^ IT)l , (7) 

Z n=— oo 

where 

Cn = I I*'* f(t)e-> nWT)t dt. (8) 

1 J-T/2 

Equation (7) may be expressed as 

h(t)=- £c n e Jlnl2 * m+aolt + l £ Cj ln{ *" T) - m]t . (9) 

A n=— oo A n=— oo 

The modulated signal expressed by (9) is composed of two doubly 
infinite spectra, one distributed about the positive carrier of radian 
frequency coo , and the other about the negative carrier of radian fre- 
quency — coo ■ It is apparent that this representation is redundant, and 
for purposes of efficient computation the number of terms can be reduced 
by a factor of two. This will be done in a way that preserves certain 
desired symmetries in the equations and allows the results to be in- 
terpreted meaningfully. First, however, it is necessary to determine the 
carrier response by introducing into (9) the expression for the trans- 
mission-frequency characteristics. The transmission-frequency char- 
acteristics given by (3) relate to baseband conditions; in order to properly 
relate this expression to carrier transmission it is necessary to translate 
these characteristics by an amount of +coo and — co . Making this trans- 
formation in (3) and substituting the results into (9) gives the following: 



git) = \ t^ C n A [» (??) + coo] a^/^V-*^' 

+ iJL^[»(r) 



T)+u ] 



(10) 

j[»(2*7r)— «0] -#[n(2ir/r)-wo] 



— I - coo I e e 

which is the expression of the response of the carrier signal to the trans- 
mission system. 

Use is now made of the fact that /(f) is an even function. The quantity 
C„ , which can, therefore, be expressed by taking twice its value over one 
half of its range, is given as, 

C„ =\f/ f(t)coan(^Jtdt. (11) 

When this, and also the fact that the attenuation must have even sym- 
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metry and the phase must have odd symmetry, are taken into account, 
(10) can be rewritten for n ^ 0, giving 



= E C n A In fe\ + cuo] cos {[«(y) + "»] l 



-7?r ) + coo 



+ 



£**[■&)- 



-77T I — WO 



COS 



[■(¥)- 



(12) 



■7F I — WO | / 

2tt n 



— COu 



Rather than computing with (12) as it stands, it will be more useful 
to extract the high-frequency carrier by the use of elementary trigono- 
metric relationships, giving 



(0 = cos u l 2 C„A \n l-Fn) + wo 



(, ° s " I y ' ~ ^ 



J + oj 



— sin coo< ^2 C n A 

n=0 



T 



WO 



> m " I y ) ; - ^ 



» + -] 



cos o>o£ 



S^[-(t 



— Wo 



2tt 



(13) 



cos-(n( — ) i - ^ | n 



(¥)-] 



+ sin w u £ YL C„A \n I — J — w 



sin <n 



r 



[•s: 



/ — ip \n I — J — o) 



Equation (13) is a general expression for the received carrier signal. 
An expression for the received response at the detector output, g(t), can 
be obtained by considering only the envelope of this carrier. This allows 
g(t) to be evaluated for any arbitrary attenuation or phase characteristic. 
In applying these results in the work to follow, the complete generality 
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of (13) will not be preserved, since a special type of transmission char- 
acteristic is being considered. By sacrificing some of the generality, how- 
ever, an advantage is gained in the facility with which the results can be 
logically organized and utilized. 

In evaluating the way a data system's performance varies as the trans- 
mission-frequency characteristics are altered, the results will be more 
lucid and systematic if the phase and attenuation are considered sep- 
arately. While these two qualities are not independent in physical trans- 
mission facilities, it is not entirely academic to treat them independently, 
since this approach is often used in equalizing facilities. 

In this study the influence of sinusoidal phase variation as given by 
(4) will be considered in some detail, and, in order to isolate the effects 
of this influence, the attenuation will be considered equal to unity over 
the frequency band of interest. (A similar study could be made of the 
effects of attenuation by specifying linear phase characteristics and 
allowing the attenuation to vary systematically.) For computational 
purposes we will use an approximate expression for r(t) by considering a 
finite number of frequency components. Applying these conditions to 
(13) and rewriting this equation for notational convenience, the carrier 
response may be expressed as 



(/) = (M + N) cos «„« + (R + S) sin a t 

= [(M + N)* + (R + Syf cos M + ¥»(*)], 



(14) 



where 

in-phase components: 



*-£c.«»p£« + 6 ■»££» + •)], 

quadrature components : 



(15) 



and 



*«> -*»"(£$$)■ («o 
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In (14), the expression for the carrier response, the envelope of the 
carrier signal that is the output of the envelope detector is given by 

y(t) - y/{M + Nf + (R + -S) 2 - (17) 

The desired expressions, which relate the output of the detector g(t) 
to the input signal /(/) (as manifest in ('„), and the three parameters b, 
t and 0, which describe the phase characteristics of the transmission 
path, are given by (15) and (17). It is these equations that form the 
basis of the computational work on the digital computer. A careful 
examination of (hem will give not only a valuable insight into the in- 
fluence of sinusoidal phase eharacteristics on AM transmission, but will 
also serve as a guide in computational work and help in interpreting 
results. It is evident from (14) that the terms M and N may be con- 
sidered as modulating a cosine carrier, and are, therefore, referred to as 
the in-phasc components of the received signal, while the terms R and S 
modulate a sine carrier and are designated as quadrature components. 
The expression for the envelope response, (17), indicates that these two 
quantities cannot be superimposed directly, but must be added at right 
angles. Also it can be seen from (15) that for b = the quadrature 
components cancel and the in-phase components add together, and the 
received signal is an exact replica of the transmitted signal. For b ?± 
but = 0, the quadrature terms cancel and the in-phase terms add to 
give a distorted replica of the transmitted signal at the receiver. For 
the case b ^ and 6 ^ 0, the received signal is a distorted replica of 
the transmitted signal that contains both in-phase and quadrature com- 
ponents. 

II can be seen from (3) and Fig. 3 that a reversal of the sign in b 
amounts to inverting the phase curve with respect to the linear phase 
line. It has been shown 1 that such a reversal in the phase curve results 
in a reversal in time of the received signal. This important fact will be 
utilized later in interpreting the performance results. Another useful 
observation that can be made from examining (15) is that a reversal in 
the sign of in these equations leaves the received signal unaltered. 

V. COMPUTED PULSE RESPONSE 

Kquations (1), (15) and (17) provide a means for computing the out- 
put of the envelope detector of an AM transmission system subject to a 
sinusoidal variation in the phase characteristics, when the periodically 
recurring raised cosine pulse of Fig. 3 is applied at the input. These 
equations were programmed for numerical evaluation on the IBM 704 
computer. In evaluating these equations, the computations were carried 
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out for a specific bit speed and specific values of b, r and 6. The results 
from these computations, however, can be given on a normalized basis. 
To do this it is necessary to express the frequency scale of the phase or 
delay curves in terms of the bit speed rather than in cycles per second, 
and also to express the time scale of the response function and the delay 
curve in terms of the bit length rather than in seconds. 

A number of normalized pulse responses are shown in Fig. 5, together 
with the normalized delay that produces the responses. The curves of 
Fig. 5(b) are for the case of symmetrical delay (6 = 0) as shown in Fig. 
5(a). For these cases, the quadrature component is zero, as would be 
expected from (15). The response curves in this figure are plots of the 
ill-phase component of the received signal. 
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Fig. 5 - (a) Normalized delay; (b) corresponding pulse response. 
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To illustrate how these normalized results can be related to specific 
cases, consider an input pulse transmitted at a bit speed of 1200 bits per 
second (bit length = 833.3 microseconds). From Fig. 5(a) the period of 
the delay curve Fo is equal to four times the bit speed or 4800 cycles, 
giving a value of r = 1/Fo = 208.33 microseconds. The curve for 6 = 8 
in Fig. 5(a) has a maximum delay excursion of br = 8 X 208.3 micro- 
seconds = 1666.3 microseconds = 2 bit lengths. The pulse response cor- 
responding to this delay is shown in Fig. 5(b) for b = 8. These same re- 
sults may also be applied to other bit speeds. For a bit speed of 2400 
bits per second (bit length = 416.7 microseconds), F = 9600, giving 
t = 104.2. Again the curve for 6 = 8 has a maximum delay of two bit 
lengths = br = 833.3 microseconds, and the same pulse response is pro- 
duced here as for the low bit speed case. 

In Fig. 6 pulse response curves are shown for conditions of asym- 
metry in the delay characteristics. This gives rise to the quadrature 
component in the received signal, and the resulting effect on the envelope 
response is shown. Here again the curves are normalized with respect to 
bit length. For the curve of Fig. 6(a), 6 = ±45°; for Fig. 6(b), = ±90°. 
These curves illustrate the fact pointed out previously: that the response 
is independent of the sign of the angle 8. 

These distorted pulse response curves indicate a degradation in 
transmission manifest in a general lowering of the pulse peak value and a 
spreading of the pulse in time. For the cases of larger distortion this re- 
sults in a considerable amount of intersymbol interference between 
adjacent pulses in a data signal. This, however, is only a qualitative 
picture of the influence of delay distortion on system performance. In 
order to make this picture more precise, it is necessary to consider a 
signal consisting of a number of pulses making up a pattern such as that 
illustrated in Fig. 2. This extension to the basic operation of computing 
pulse response was achieved by combining the received response of a 
single pulse (before envelope detection) in accordance with the desired 
pattern. In doing this, the in-phase components from adjacent pulses 
were added together. Similarly, the quadrature components were added 
together, and the resulting envelope response was determined in ac- 
cordance with 

g(t) = V[S(M + N)Y + [2(B + S))*, (18) 

where M, N, R and S are defined by (15). By this method, the results 
previously discussed can now be extended so that the output response 
for pulse patterns can be obtained. 

An example of pulse-pattern response showing the resultant in-phase 
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Fig. 6 - Pulse response curves. 

and quadrature components and the envelope response is given in Fig. 7. 
In this figure a pattern of binary characters is shown together with the 
corresponding output response curve for conditions b = 1, br = 0.75 bit 
length and 6 = 90°. In order to recover the binary information from the 
signal, the envelope response is examined at the sampling-time points to 
determine if the value of the curve at each of these points is less than or 
greater than the value of the slice level. The distance from the signal to 
the slice level for a given time point (as indicated by the arrows) is a 
measure of the noise margin for that point. The amount by which the 
noise margin is reduced by the transmission characteristic can be de- 
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Fig. 7 - Basic method for evaluating performance. 

termined by comparing these values with the corresponding values for 
an undistorted or reference output resulting from ideal transmission 
characteristics. 

This procedure makes possible a quantitative evaluation of data 
system performance. It is evident, however, that for longer pulse pat- 
terns this procedure would soon become tedious and impractical if car- 
ried out in the manner shown. In order to alleviate this difficulty, the 
basic method for evaluating performance, as illustrated in Fig. 7, was 
included in the computer operation for any desired pulse pattern up to a 
certain maximum length. As the delay distortion increases, the range of 
intcrsymbol interference of a single pulse becomes larger; it was found 
necessary in this analysis to include cases for which the intersymbol 
interference extended over as many as nine bit lengths. In order to take 
this influence into account it was necessary to use a pattern that con- 
tained all possible combinations of binary characters that could occur in 
a nine-bit word. This results in a maximum pattern length of 2" = 512 
bits. A binary pattern of this type was used in the computer simulation 
work that was carried out. 



VI. EYE PATTERN' 

One way of portraying this procedure for evaluating data system 
performance for a long sequence of pulses is by means of the "eye" pat- 
tern. The manner in which this pattern is made can be illustrated by 
dividing the output response signal from the system into a number of 
segments each of duration equal to two bit lengths and plotting all these 
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segments over the same two-bit time intervals. Such a plot for a reference 
or undistorted eye is shown in Fig. 8(a). Here the eye pattern is made 
up of traces that are straight lines at the top and bottom of the eye, 
corresponding to a sequence of ones and zeros, respectively, and traces 
that have sinusoidal shape in between, representing transitions from one 
to zero and zero to one. Since the pulses are undistorted, all traces of a 
particular type fall on top of one another, and the resulting eye pattern 
is made up of well-defined lines. In this case, the opening or aperture A, 
of the eye is equal to A , which is the maximum value of a single pulse. 
The quantity A is a measure of the noise margin and will be used as a 
criterion for data system performance. 

For conditions of transmission distortion, the resulting eye pattern is 
illustrated by Fig. 8(b). Here the various individual traces have been 
distorted and no longer form well-defined lines of transition, but are 
spread out into bands. This results in a reduction of the aperture, which 
can be expressed numerically as an impairment in the system perform- 
ance : 



transmission impairment = 20 logio -7- • 

A 



(19) 



There are two factors that determine the validity with which such an 



(a) 




(b) 



Fig. 8 - (a) Reference eye pattern; (b) received eye pattern with transmission 
distortion. 
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expression of performance describes an actual system. First, it must be 
acknowledged that the effective opening of the eye may change from 
one trace to another, and the quantity A as shown in Fig. 8(b) is a mini- 
mum bound of this opening. Since an opening of this magnitude will 
occur relatively infrequently this quantity is then a pessimistic estimate 
of the system performance; getting an exact measure would involve 
taking account of the statistical distribution of the quantity A. The 
probability of an error occurring would then be determined by con- 
volving this distribution with the distribution of noise. Such a refine- 
ment is rather complicated, and may not be necessary except in cases of 
very small aperture, since, for the larger values of A, the traces across 
the eye are rather densely packed. 

The other fact to be considered is that in this simulation the degrada- 
tion in performance due to timing recovery is not included. These results 
then specify the performance under conditions of ideal timing recovery, 
or for systems operating with ideal transmitted timing. Again, this is a 
good approximation to system operation with timing recovery, except 
for small aperture values, where timing jitter is unavoidable. 

VII. EXPERIMENTAL VERIFICATION OF COMPUTER SIMULATION 

It is of interest to determine how closely the performance results from 
this simulation procedure check with measured performance from an 
actual data system. The aperture value of an AM double-sideband data 
system operating in the voice band frequency range was measured at 
several bit speeds over a given transmission channel. The transmission 
characteristics of the channel were then applied to the computer simula- 
tion of the data system for the same bit speeds. The results from these 
two tests expressed in terms of aperture of the eye pattern are shown in 
Fig. 9. The fact that the computer simulation performance is slightly 
better over the range of bit speeds tested is to be expected, because of 
the assumption of idealization in the modulator and detector in making 
the simulation. These results tend to confirm the validity of the simula- 
tion technique as a means for evaluating data system performance. 

VIII. RESULTS 

A criterion for data system performance has been specified as A , the 
aperture of the eye pattern formed by the received signal of a data 
system. The results of computing performance based on the data system 
simulation are presented here as a portrayal of the manner in which A 
varies with a systematic variation of the parameters describing the 
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Fig. 9 - Comparison of computer simulation performance with experimental 
verification. 

transmission characteristics of the system. To make these results more 
general, a normalized parameter 5, will be defined as 8 = t X bit speed. 
Results showing the manner in which the quantity .4 varies with b, 8 
and 8 arc given in Figs. 10 through 20, where b is expressed in radians 
and 6 in degrees and 8 is a quantity equal to the reciprocal of the period 
of the sinusoidal phase characteristics normalized with respect to bit 
speed. 

In these figures, performance curves are plotted with the dependent 
variable b ranging over values sufficient to cause the aperture to vary 
from a reference value of unity to complete closure. The quantity 5 
varies from 0.2 bit length to 1 .4 bit lengths. The degree of asymmetry 
is specified by the angle 0, which varies as a parameter from zero to 
90 degrees. It is an apparent feature of each set of curves that the aper- 
ture value is reduced more rapidly, as b increases, for asymmetric phase 
conditions than for symmetrical conditions. This is an evidence of the 
influence of the quadrature component on data system performance. 
From these curves, it can also be seen that the relative influence of the 
quadrature component is much more pronounced for small values of 8 
than for larger values. 

While the performance curves cover only one quadrant of the sinusoi- 
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Fig. 10 - Plot of aperture vs. 6 for 5 = 0.2 = t X bit speed. 
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Fig. 11 - Plot of aperture vs. b for 5 = 0.3 = t X bit speed. 
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Fig. 12 - Plot of aperture vs. b for 5 = 0.4 = t X bit speed. 
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Fig. 13 - Plot of aperture vs. b for 5 = 0.5 = t X bit speed. 
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Fig. 14- Plot of aperture vs. b for 5 = 0.6 = t X bit speed. 
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Fig. 15- Plot of aperture vs. b for 6 = 0.7 = r X bit speed. 
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Fig. 16- Plot of aperture vs. b for 5 = 0.8 = t X bit speed. 
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Fig. 17 - Plot of aperture vs. b for 8 = 0.9 = t X bit speed. 
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Fig. IS- Plot of aperture vs. b for 5 = 1.0 = t X bit speed. 
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Fig. 19- Plot of aperture vs. b for 5 = 1.2 = r X bit speed. 
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Fig. 20 - Plot of aperture vs. b for 8 = 1.4 = r X bit speed. 

dal phase curve, the other three quadrants are duplicates of the curves 
shown. This can be seen from the facts previously stated, which indicated 
that a change in the sign of 8 leaves the output response unchanged, and 
that a change in the sign of b reverses the time variable of the response 
and therefore reverses time in the eye pattern but does not change the 
aperture value. 

The information given on the performance curves of Figs. 10 through 
20 provides a means for evaluating the aperture of a data system for a 
wide range of sinusoidal phase characteristics, either by direct means or 
by interpolation. These results can be displayed in another manner, 
which gives a more concise picture of how the performance varies. To do 
this the performance information for a particular value of 8 is plotted in 
polar form with b and 6 as the radial and angular coordinates respec- 
tively. Such a plot is shown in Fig. 21 for 8 = 0.75, with contours of 
constant aperture shown as a parameter. The outer contour has a value 
of A = and is the locus of points for which the aperture is completely 
closed. Any sinusoidal phase characteristic whose period corresponds to 
this value of 8 can be located on the polar diagram of Fig. 21 as a point. 
If the magnitudes of b and are such that the point lies within the outer 
contour, the aperture value can be determined either directly or by inter- 
polation. It is now possible to plot a path on this diagram showing how 
the value of A varies for changes in the values of b and 9. 
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This type of plot can be made for other values of 8. By combining such 
plots, the information given in Figs. 10 through 20 results in a very use- 
ful three-dimensional model with coordinates b, 8 and 0, which depicts 
the performance of the data system. Such a model is shown diagram- 
matically in Fig. 22. In this three-dimensional space any sinusoidal 
phase (or delay) curve is represented as a point. A number of curves 
that illustrate the particular delay that occurs at various points are shown 
by Fig. 22(a), (b), (c) and (d). The surface of the model indicated by the 
dashed line in Fig. 22 represents the locus of points for which the aperture 




Fig. 21 - Polar diagram of b vs. for 8 = 0.75 with aperture as a parameter. 
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value is zero. For points inside this surface the aperture is greater than 
zero. Other surfaces within this model can be described by connecting 
all points having a common value of A. The points lying along the verti- 
cal axis represent conditions of zero distortion, or maximum aperture 
value A. A photograph of a performance model constructed from the 
numerical results obtained from the computer simulation is shown in 
Fig. 23. 



IX. APPLICATION OF RESULTS TO ACTUAL TRANSMISSION FACILITIES 

In the simulation of the data system that has been presented the phase 
characteristics have been assumed to be sinusoidal in shape, for reasons 
of mathematical convenience. Some rather general results have been 
given in terms of the parameters that describe the phase characteristic. 
In this section several specific transmission facilities will be considered in 
the light of this simulation technique. 

The validity with which the sinusoidal phase (or cosine-shaped delay) 
represents actual transmission characteristics will be investigated, and 
computed performance resulting from data system operation over these 
facilities is given. 

In Fig. 24, the delay characteristics of four message facilities are shown 
by the solid lines, together with a cosine approximation to these delay 
curves given by the dashed lines. It can be seen from Fig. 24 that the 
cosine-shaped delay curves give a very good approximation to the actual 
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Fig. 22 - Three-dimensional model of data system performance. 
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Fig. 23 - Performance model constructed from numerical results of computer 
simulation. 

delay over the frequency range of 1000 to 2000 cycles per second, the 
maximum departure of the approximating curve from the actual facili- 
ties being less than 100 microseconds over this frequency range. The 
approximating curves can then be specified by the three parameters b, 
t and (or when the bit speed of the data system is given by b, 8 and 0). 
These values are given in Table I, for the cases of a data system operating 
with an 1800-cycle carrier, and at bit speeds of 1200 and 1500 bits per 
second. The resulting aperture values and transmission impairment fig- 
ures produced are also given in Table I. Aperture values were determined 
by direct computer computation; these values can also be determined 
to a good approximation by interpolation from Figs. 10 and 11. 

The method shown here for determining the impairment values in 
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Fig. 24 - Delay characteristics of four message facilities. 



Table I illustrates how data system performance can be established for 
actual transmission facilities from the results of the computer simula- 
tion. This procedure can also be used to determine how performance 
varies with some of the parameters that are important in the design and 
engineering of data systems, such as bit speed, location of carrier fre- 
quency and length and type of transmission facility. For example, as 
the bit speed is increased from 1200 to 1500 bits per second, the im- 
pairments of the facilities are increased by the amount shown. (The 
loaded cable changes from marginal operation to complete failure.) 

The poor performance of the loaded cable facility is due to the fact 
that the location of the carrier frequency results in an asymmetrical 
disposition of the delay curve about the carrier, resulting in a consider- 
able amount of quadrature component. If the location of the carrier were 
changed from 1800 to 1500 cycles per second, the angle 6 would be 
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Table 


I 














ra- 
dians 


T, 

H sec 


o, 

de- 
grees 


1200 bits per second 


1500 bits per second 


Facility 


5. 

bit 

lengths 


A* 


Trans- 
mission 
impair- 
ment, 
db 


bit 
lengths 


A* 


Trans- 
mission 
impair- 
ment, 
db 


400-mile H44 loaded cable 

4 links of N carrier 

5 links of K carrier 
8-repeater sections of 

V.F. open wire 


9.26 
11.8 

9.89 
11.2 


194 
163 
190 
161 


43.8 

15.8 

-3.1 

3.6 


0.233 
0.196 
0.228 
0.193 


0.092 
0.60 
0.67 
0.76 


20.7 
4.4 
3.5 
2.4 


0.292 
0.244 
0.285 
0.242 


0.38 
0.50 
0.61 


8.4 
6.0 
4.7 



* These aperture values were determined by direct computation. They can also 
be obtained by interpolation from Fig. 10 and Fig. 11, with the exception of the 
loaded cable (which has an aperture value of either zero or very nearly zero). 

changed from 43.8° to 22.8°. For a bit speed of 1200 bits per second, this 
gives an aperture value of 0.41 (by interpolation from Figs. 10 and 11) 
and a resulting transmission impairment of 7.8 db. 

Performance for many other transmission conditions of interest can 
be determined with similar ease by this same procedure. 



X. CONCLUSIONS 

It is apparent that the data system simulation technique that has 
been developed and used here for evaluating performance is a powerful 
and valuable tool. Wherever possible, experimental and theoretical veri- 
fications of the results obtained from this simulation have been made. 
These tests indicate that this method is valid and accurate. The per- 
formance curves given have been used to illustrate the ease and effec- 
tiveness with which many of the complex problems relating to delay dis- 
tortion can be answered in a rather general way. 

Several uses can be made of this simulation technique. First, it pro- 
vides a means of obtaining quantitative answers to a broad category of 
problems relating to data transmission, such as the influence on per- 
formance of the variation of bit speed, the location of carrier frequency 
and the length and type of transmission facility used. In addition, a 
performance model based on a sinusoidal phase characteristic and nor- 
malized with respect to bit speed has been produced for the data system 
under consideration. The surface of this model describes the boundary 
(in terms of variation of phase distortion) within which the system can 
operate. Similar models can be established for other types of data sys- 
tems to provide a means for comparing the performance on a much 
broader basis than is presently possible. 
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The investigation presented here was directed to the problem of deter- 
mining what influence phase (or delay) distortion has on data system 
performance. It is clear that the simulation method developed could be 
extended to include other problems of data transmission, such as a 
general study of the influence of the attenuation characteristics, the 
problem of recovering timing and the influence of noise on system trans- 
mission. Additional work in this general area is now being carried out. 
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